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The Schwarzschild black hole can be viewed as the special case of the marginally bound Lemaitre- 
Tolman-Bondi models of dust collapse which corresponds to a constant mass function. We have 
presented a midi-superspace quantization of this model for an arbitrary mass-function in a sepa- 
rate publication. In this communication we show that our solution leads both to Bekenstein's area 
spectrum for black holes as well as to the black hole entropy, which, in this context, is naturally 
interpreted as the loss of information of the original matter distribution within the collapsing dust 
cloud. 
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I. INTRODUCTION 

Due to the singularity theorems of Penrose, Geroch, 
and Hawking pj and the many examples that have been 
found in various models of gravitational collapse it 
has now come to be generally accepted that the gravi- 
tational collapse of a sufficiently massive star will lead 
eventually to the formation of a black hole. 

In one particular class of models, the LeMaitre- 
Tolman-Bondi (LTB) models of spherical, inhomo- 
geneous dust collapse, the picture that emerges is sim- 
ple and general enough to be interesting. These mod- 
els are all represented by a solution of Einstein's equa- 
tions with pressureless dust described by the stress ten- 
sor T^ v = eU^Uv, where e(t,r) is the energy density of 
the collapsing dust cloud. The solution is characterized 
by two arbitrary functions of the coordinate r, the mass 
function, F(r), and the energy function, /(r), the former 
being related to the initial energy density distribution in 
the collapsing cloud and the latter to its initial veloc- 
ity profile. The collapsing dust cloud is imagined to be 
made up of shells, labeled by r, which successively col- 
lide at late times into a central singularity. The classical 
solution is given by 
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where R is the physical radius, R' represents a derivative 
w.r.t. r and R* represents a derivative w.r.t. the dust 
proper time, r. 

An important subclass of these solutions, in which the 
energy function is taken to be exactly vanishing, de- 
scribes the marginally bound models. In this case, using 
a freedom to scale the initial distribution, the general 



solution of (1.1) can be cast into the form 



R*(r,r) = r* - -VF^jr. 



(1.2) 



The epoch R = describes a physical singularity. Suc- 
cessive shells, labeled by r, collapse into this singularity 
at proper times given by 



r(r) 



2ri 



(1.3) 



Various models, with remarkably different behaviors, 
are obtained from different choices of the mass func- 
tion, F(r). Some models lead to the formation of a 
naked singularity for certain initial data and of a black 
hole for other initial data. One particular model, with 
F(r) — 2M (a constant), does not describe dynamical 
collapse at all but rather a static Schwarzschild black 
hole of mass M . This is most easily seen by performing 
the following trans form ation from the Tolman-Bondi co- 
ordinates (t, r) in (1.1) to curvature coordinates, (T,R), 
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to recover the usual static, Schwarzschild form of the 
metric. T is the Schwarzschild (Killing) time. 

Black holes are extremely interesting objects. Al- 
though their temperature is exactly zero Kelvin in clas- 
sical general relativity, Bekenstein proposed 0] that they 
have a temperature and an entropy and should be treated 
as thermodynamic systems. Their temperature and en- 
tropy are known from semi-classical arguments |Bj to 
be fundamentally quantum mechanical, yet the precise 
quantum origin of these properties remains shrouded in 
mystery despite the many interesting proposals || that 
have been made in recent years. Ultimately, black holes 
are the classical end states of collapse and it should be of 
considerable interest to understand these properties from 
a bona fide microcanonical ensemble of quantum states 
constructed from a collapsing matter distribution. In this 
communication we will examine the quantum mechanics 



of black holes from the vantage point of a midi-superspace 
quantization of the LTB models described above. At least 
two problems are eliminated by this approach: (i) the 
proper time of the collapsing dust enters naturally in the 
dynamical equation, eliminating the so-called "problem 
of time" that usually plagues the canonical quantization 
of gravity and (ii) there are no operator ordering ambigu- 
ities, the configuration space is flat and an inner product 
may be defined unambiguously on the Hilbert space so 
that the usual probabilistic interpretation of quantum 
mechanics is retained. 



II. QUANTIZATION OF DUST COLLAPSE 

In a separate publication M we have presented the 
canonical quantization of the marginally bound models 
for a variety of mass functions, including those that admit 
both black holes and naked singularities. After a series 
of canonical transformations, performed in the spirit of 
Kuchaf 's remarkable reduction of static spherical geome- 
tries pi, it becomes possible to describe the phase space 
by the dust proper time, r(r), the physical radius, R(r), 
the mass function, F(r) and their canonical momenta, 
P T (r), Pji(r) and Pf{t). The collapse is reduced to two 
classical constraints, one of which, the momentum con- 
straint, 



t'P t + R'P R + F'P F 



0. 



insures spatial diffeomorphism invariancc. 
the Hamiltonian constraint, 

(P T + F'/2) 2 +TP 2 R -^^0, 



(2.1) 
The second, 



(2.2) 



is responsible for (p roper) time evolution (we have set 
T = 1 — F/R in ( |2.2[ )). The mass function, F(r), defines 
the particular collapse model being considered. This is 
best seen in the Lagrangian description in which the gen- 
eral s olu tions of Einstein's equations are given by (1.1) 
and (1.2). These solutions show that F(r) is related to 
the proper energy of the collapsing cloud. (F'(r) is the 
energy per unit coordinate cell, d 3 x, and is requ ired to 
be positive definite.) As a consequence of ( |l . l| ) , F(r) 
determines e(r, r) and vice- versa, apart from a choice of 
scaling, e(r, r) determines F(r) according to 



F(r) 



e(r',0)r' dr', 



(2.3) 



where we have used the scaling implicit in (L2). Both 
F{r) and e(r, r) are externally given functions that de- 
fine the collapse. This comes about because the energy 
density, e(r, r), is a Lagrange multiplier, as seen from the 
form of the dust action, 



~ J d 4 x^ e(x) [g^U 13 + 1] 



(2.4) 



that reproduces Einstein's equations. As a Lagrange mul- 
tiplier, e(r, r) enforces "time-like dust" i.e., it requires 
the dust world-lines to be time-like geodesies. Thus F(r), 



like e(r, r), can be chosen arbitrarily, but a choice spec- 
ifies an initial energy density distribution and thus a 
collapse model. In the quantum theory, F(r) acts as 
a weight in the potential term, which identifies the col- 
lapse model being quantized. Each model corresponds to 
stellar collapse under certain initial conditions. 

The DeWitt supermetric, considered on the effective 
configuration superspace, ( r, R ) , is non-degenerate, can 
be read off directly from ( |2.2| ) and is found to be flat, 
positive definite when T > and indefinite when T < 0, 
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(2.5) 



It happens that, for the spherical geometries being con- 
sidered, T — or R — F describes the apparent horizon 
while, for the black hole geometry, it is the event hori- 
zon. When T ^ 0, the supermetric can be brought to a 
manifestly flat form by the coordinate transformation 



R* = ± 




(2.6) 



The next step is to turn the canonical momenta into co- 
variant functional differential operators (covariant w.r.t. 
the supermetric, jab) 0, according to 
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(2.7) 



which act on a state functional, '5 [r, R, F] . Then defining 

*[t,R,F] =e~^r F '^ T ^ dr ^[ T ,R,F}, (2.8) 

the functional \P is seen to obey the (Wheeler-DeWitt) 
equation, 
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(2.9) 



which is similar in form to a Klein-Gordon equation for a 
scalar field with a potential. Because the metric in ( 2~5| ) 
is positive definite in the region R > F, corresponding, 
in the collapse geometries being considered, to the region 
outside the horizon, and indefinite when R < F, the func- 
tional equation is elliptic in the former and hyperbolic 
latter regions. It is convenient to write this equation as 
a functional Schroedinger equation in (r, i?» ) by taking 
the square-root of (2J3) as follows, 
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*[t,R,F], 



(2.10) 



where the negative sign within the square root refers to 
the region R > F and the positive sign to the region 
R < F. Thus the dust proper time may be identified as 
the time variable as was originally proposed by Kuchaf, 
Torre and Brown pjQ|, H ]. 

Any solution of (Z.lOjrnust obey the diffeomorphism 
constraint 



5t + *5R* + SF 



$[r,R,F}=0 (2.11) 
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and we will first provide an ansatz for a solution to this 
equation: as long as F'(r) is not everywhere vanishing, 
take 



V[t,R,F] = exp 



drF' (r)W(T(r), R(r), F(r)) 



(2.12) 



where W is an arbitrary, complex valued function of its 
arguments (and not t heir derivatives) which is to be de- 
termined from ( [2.10 ). To see that this is indeed a so- 
lution, we note that the integrand in the exponent is 
clearly a spatial density because, while r(r), R(r) and 
F(r) are spatial scalars, F' (r) i s a density. It follows 
that the wave- functional i n (2.12 ) will obey the momen- 
tum constraint. Indeed, (2.11) simply requires th at W 
admit no explicit dependence on r. Our ansatz (2.12) 



is not unique, of course, but it is guided by a physical 
consideration, namely that in the classical theory F'(r) 
is related to the proper energy density of the collapsing 
shells of dust. 



Together, equations ( 2.1C ) and (2.11) define the quan- 
tum theory whose inner product is given by the func- 
tional integral, 

(*i,*a) = JvR^\^ 2 = J Y[dR*(r)*l*2. (2.13) 



As pointed out in fujl , the inner product defined in this 
manner ensures the hermiticity of the momentum, P*, 
conjugate to R„. The norm of a quantum state under 
this scalar product is formally r independent provided 
t hat h defined in ( 2.10 ) is self- adjoint. It is clear from 
( |2.9| ) that this operator represents the proper energy of 
the shell labeled by r, the total energy of the system 
being simply 



H 



drh(r). 



(2.14) 



It is sclf-adjoint only in the linear sub-space in which the 
operator =fP? + F' 2 /4!F admits positive eigenvalues. 



III. STATIC BLACK HOLES: BOUND STATES 
AND MASS QUANTIZATION OF A SINGLE 
SHELL 

We now apply the above considerations to the static 
black hole. As we have argued, in the quantum theory 
as in the classical theory, the mass function is externally 
specified and defines the collapse model being considered. 
In order to specify it in a way that is most consistent 
with the classical solution (F — 2M), we consider a mass 
function describing a single spherical shell of total mass 
M . Spherical shells have been used often in studies of 
gravitational collapse and cosmology jl2j . The shell label 
is immaterial, but for concreteness let us suppose it is 
r = 0. Then we define the shell by the mass function 



F(r) = 2M9(r), 



(3.1) 



where 9(r) is the usual step function. Therefore, 

F'(r) = 2M5(r), (3.2) 

which vanishes everywhere except on the shell itself 
where it contr ibutes infinitely to the shell self-energy. 
Equation ( 2.12Q tells us that the problem of a single shell 
is essentially quantum mechanical, 



*[r, R, F] = exp [MW(r, R, M)\ 



(3.3) 



where r = t(0),R = P(0) and F(0) = 2M represent, 
respectively, the proper time, the radial coordinate and 
the total mass of the single shell. The Schwarzschild 
black hole as a single shell has but o ne degree of freedom 
and the functional equation ( |2.10 ) tur ns i nto an ordi- 
nary Schroedinger equation. F'(r) in (3.2) contributes 
5(0) 2 to the potential term in ( |2.10|) . We take this to be 
vanishing, in keeping with DeWitt's regularization [fl3|| . 
This subtracts the shell's infinite self-energy. It should 
be remarked, however, that this regularization scheme is 
not unique and that there have been attem pts to define 
a quantum theory in which 5(0) ^= (sec fl4||). There 
is, unfortunately, no general consensus on how the coin- 
cidence limits in the Wheeler-DeWitt equation are to be 
treated and for the present we assume DeWitt's regular- 
ization. We then have 



<9* 



(3.4) 



As is normal in the quantum theory, we as sume that the 
wave function is We take P* in (2J3) to range over 

the entire real line extending, in the interior (R < F), 
from -irM to +7rM, 



P* =± 



-y/R(2M-R) 

R-M 



M tan" 
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(3.5) 



(its magnitude is the radius of the wormholc throat) and, 
in the exterior [R > F), from ±7rM to ±oo, 



P* 



± 



y/R{R - 2M) 



M\n[R - M + y/R(R-2M)} 
MlnM + irM] . 



(3.6) 



The classical singularity occurs at P* = 0. 

The stationary states of the black hole are now quite 
easily described: in the interior they are a superposition 
of ingoing and outgoing waves, 

*[r,i?»] =A ± e- lM(T±R '^ (3.7) 
and in the exterior, they are exponentially decaying 

*[t, R»] = Be- iM{T - lR '\ R* > nM 

= Ce - lM{T+lR '\ P*<-ttM. (3.8) 

Matching the wave function and its derivative at the hori- 
zon, one finds that the energy (mass) squared of the black 
hole is quantized in half integer units, 
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(3.9) 
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where M p is the Planck mass. This is the Bekenstein 
area spectrum 11. A similar result was reported earlier 
by us in |l5| (although our construction in that work 
was not aided by the collapse model we are considering 
here) with the difference that the energy of the ground 
state was found to be exactly zero. On the contrary, the 
minimum shell mass is here found to be M p / y/2. The rea- 
son for this discrepancy is that the wave-function in fl5| l 
was taken to be identically vanishing outside the horizon. 
This condition was too strong and in fact unnecessary. It 
will be seen that the non- vanishing ground state energy is 
required for the interpretation of the black hole entropy 
that follows. 



It is worth seeing how the choice, (3T), reproduces the 
constraints for massive dust shells that have been used by 
various authors in the past ]l6}| , |l7j| . These constraints 
were expressed in the phase space constructed from the 
metric coefficients in 

ds 2 = N 2 dt 2 - L 2 {dr - N r dt) 2 - R 2 dfl 2 (3.10) 

and their conjugate momenta, and were derived from a 
square-root form of the matter action. They take the 
form, 



H d 



H r . 



V- 



pS(r — f). 



(3.11) 



In this system, (t, L, R, P t , Pl, Pr), the dust constraints 
obtained from the quadratic action ( |2.4D , upon eliminat- 
ing the energy density, are (see, for example, or fll|| ) 



H" 



u 2 

^ r 

L 2 



K = -U r P Tl 



(3.12) 



w here U r = —t'. It is easy to see that the constraints 

f. 12 ) turn into (3.11) by taking P T = mS(r — f) — 
'(r)/2, as we have above (f replaces "0" as the shell 
label). The dust supermomentum is = —P T U r = 
mr'<5(r — f ) = —p5(r — f), where we have set p = mU r . 
The gravitational part of the constraints, of course, has 
the usual form. 



IV. STATIC BLACK HOLES: THE 
SUPERPOSITION OF SHELLS AND THE 
STATISTICAL ENTROPY 

Black holes are not generally formed by the collapse of 
a single, infinitesimally thin, shell of matter but by many 
shells, progressively colliding into a central singularity. 
The static black hole must be viewed as the final state of 
their collapse, i.e., when all shells have collided with the 
central singularity (see figure 1). 




constant 
horizon 



all shells 
"inside" 



Fig. 1: Collapsing shells form a black hole 

Let us therefore consider the description of N shells. We 
tak e a simple minded generalization of the mass function 
in (0), 



3=1 
N 



F'{r) = 2Y,»Ar-r 3 ). 

3=1 



(4.1) 



Inserting F'(r) into ( 2.12|) s hows that the wave- functional 
is now a product state ||l8| , 



N 



(4.2) 



fe=i 



over N wave-funct ion s, one f or ea ch shell. The coordinate 
is given by (3.5) and ( |3.6|) and is determined not 
simply by the mass, Hk , of the shell (as it was for a single 
shell) or the total mass, M , of all the shells, but by the 
total mass contained within it, 

k rr k 

2M k = 2 V fij = / drF'{r) = F{r k ). (4.3) 

Information about the gravitational interaction between 
shells is thus seen to be encoded in i?*fc, which is 



+ Mfetan 1 
inside, and 



-y/R k (2M k -R k ) 

Rk - M k 



y/R k (2M k -R k ) 



7rM fc 



(4.4) 



R,„ 



± 



y/R k (R k -2M k ) 



M k ]n[R k - M k + y/R k (R k -2M k )} 
M k \nM k +TTM k ] 



(4.5) 



outside. We see, not surprisingly, that the collapse of 
each shell is sensitive to the precise mass distribution 
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among the shells. Our system of N shells has N de- 
grees of freedom and each wave-function obeys th e or - 
dinary one dimensional Schro edin ger e quation in (3.4), 
with solutions given by (|3.7|) and (3.S), again eliminat- 
ing the infinite shell self-interaction in the potential. A 
straightforward application of the boundary conditions 
appropriate to each shell gives the following quantization 
condition for the states of shell k, 



HkM k 



(4.6) 



S = fcln 



{N/2 + (M/M p ) 2 - 1)! 
(N- 1)\[{M/M p f -N/2]\ 



(4.10) 



This result depends on the number of shells that have col- 
lapsed to form the black hole. However we do not know 
the number of shells that formed the black hole because 
the macrostate is defined only by the black hole mass, 
so this number should be independently determined. For 
an eternal black hole in equilibrium it is natural to max- 
imize the entropy with respect to N. When both TV and 
M/Mp are large, one readily finds 



These conditions, if applied recursively, show that the 
mass of shell k is determined by k quantum numbers. 
Thus the total mass depends on N quantum numbers for 
a black hole formed out of N quantum shells. 

We are now able to understand the origin of the black 
hole entropy in tangible terms. The appearance of N 
quantum numbers means a quantum black hole is not 
simply described by its total mass. Such a description ig- 
nores the manner in which the mass is distributed among 
the shells. The entropy counts the number of distribu- 
tions for a given total mass. 

The reasoning leading up to ( |4.6|) cannot be completely 
extended to the static black hole region. This is because 
it relies on a certain classical ordering of the shells, i.e., 
shell one is "inside" shell two, which is "inside" shell 
three, etc. But this ordering makes sense only as long as 
the shells have not collided into the central singularity 
(R = 0). Because, in the static region all the shells have 
collapsed to the same physical point, all information of 
the original spatial mass distribution is completely lost 
and we know only that there is one horizon (of physical 
radius R = 2M) for all the shells. Consequently, the 



mass condition in (4.6) should be replaced by the simpler 
relation 



fi k M 



nk 



I) 



(4.7) 



where M is the total mass of the hole. Again, the total 
mass (squared) of the hole continues quantized, 



M 



N 1 

2> = mE 



- I \f~ 

2 / v> 



(4.8) 



but now in integer as well as half integer units. 

The problem of counting the number of distributions 
will be recognized as one from elementary statistical me- 
chanics texts. It is precisely the problem of N simple 
harmonic oscillators into which one wishes to distribute 
a total number, say Q, of quanta. (Alternatively, ask for 
the number of ways in which N integers may be added 
to give another integer, Q.) If we knew the number of 
shells that went into the black hole's making, the answer 
would be 



n = 



(N + Q-l)\ 
(N-1)\Q\ 



(4.9) 



In this case we have Q = (M/M p ) 2 - N/2, from (| 
The statistical entropy would then be (exactly) 



TV w — ( — 



giving, to leading order, 



S ~ k In p 



M 



(4.11) 



(4.12) 



where 



( 1 + ts) <1+ * ) 



_2_\ n/5 (-. l_\V VE) 

V5 I 1 V5j 



2.618. 



(4.13) 



If we now insert (4.11) into ( |4.8| ) we find that the to- 
tal mass of the black hole is quantized according to the 
relation M 2 m uM*/(l - 1/V$) for v £ N, so that the 
entropy can be given as 5* « vkhxp/{\— 1/ y5) . The aver- 
age mass square of each quantum shell is y/5M p /2 <~ M 2 , 
therefore each shell may be thought of as contributing 
about one Planck mass to the total mass of the black 
hole. In other words the shells, on average, are all virtu- 
ally in their ground states. 



V. CONCLUSIONS 

We have shown how the black hole area quantization 
and entropy can be understood in terms of the collapse 
of shells of matter. The entropy encodes the loss of in- 
formation of the matter distribution among the shells in 
the final state. Shell states do not in general vanish at 
the classical singularity, so the classical singularity is not 
prohibited by the quantum theory. However, it is inter- 
esting that the wave function at the singularity either 
vanishes or has vanishing derivative. 

What are the solutions when the matter distribution 
is continuous? Various continuum models lead to the 
formation of both naked singularities as well as black 
holes therefore the quantization of the continuum 
models may be expected to shed some light on Hawk- 
ing radiation and the Cosmic Censorship Hypothesis. A 
solution for a general differentiable mass function was 
given in Q). However, a realistic collapse would involve 
one or more regions of differing mass functions and while 
the mass function must be continuous across the bound- 
ary between the regions, it may not necessarily be dif- 
ferentiable. The wave-functional must be appropriately 



5 



matched, i.e., required to be both continuous and differ- 
cntiablc at every such boundary. Specific models will be 
described elsewhere. 
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